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Wave propagation through waveguides, quantum wires or films with a modest amount of disorder 
is in the semi-ballistic regime when in the transversal direction(s) almost no scattering occurs, while 
\^ ' in the long direction(s) there is so much scattering that the transport is diffusive. 

^\ , For such systems randomness is modelled by an inhomogeneous density of point-like scatterers. 

0\ ■ These are first considered in the second order Born approximation and then beyond that approxima- 

tion. In the latter case it is found that attractive point scatterers in a cavity always have geometric 
f— I , resonances, even for Schrodinger wave scattering. 

' In the long sample limit the transport equation is solved analytically. Various geometries are 

considered; waveguides, films, and tunneling geometries such as Fabry-Perot interferometers and 
double barrier quantum wells. The predictions are compared with new and existing numerical data 
and with experiment. The agreement is quite satisfactory. 

pacs 

71.55 Jv (Disordered structures) 
73.20 Dx (Electrons in low dimens. struct.) 
73.50 -h (Electron transport in thin films) 
42.81 Dp (Propagation in fiber optics) 
73.40 Cg (Contact resistance) 
^) ' 73.40 G (resonance tunneling) 

S I As is empirically known from the ancient development of music instruments, cavity resonances determine the 
■ transmission of waves through devices which have dimensions of the order of the wavelength. The resonances give 
(~| ' rise to transmission peaks of diverse systems such as flutes, organ pipes, Fabry-Perot interferometers, electronic 
O nanostructures [Q, and electronic waveguides. In these systems the transmission of waves can drastically increase if 
^ ' the wavevector of the incoming waves allows for a new mode to be resonant. In the case of a very pure or very small 
^ , cavity, impurity scattering can be neglected and the transmission is said to be ballistic. Ballistic transport has been 
shown to occur in various systems, including quantum point contacts and narrow optical slits 

For wave propagation through waveguides or quantum wires with a modest amount of disorder, several regimes 
occur. For rather clean systems one still has ballistic transport of essentially unscattered waves. In the limit of dirty, 
■ ■ ' but non- absorptive and phase coherent systems, the intensity mainly diffuses through the system. As pointed out 
recently by one of us, for long wires or thin films there is a third regime, the semi-ballistic regime. Q Here disorder 
is large enough to cause diffusion in the long direction(s), but small enough to maintain ballistic motion in the narrow 
directions. In the present work we shall focus on this regime. 

For electronic systems the conductance can be expressed by the Landauer formula, 

_ 2e^ \ ^ ^flux 

ab 

in terms of the flux transmission coefficients T^^^ of the system, where a and b stand for the incoming and outgoing 
channels, respectively. In the ballistic regime there is no channel-to-channel scattering and the transmission coefficients 
are diagonal in mode space. They are of order unity for the propagating modes (the low order cavity resonances that 
can be excited at the energy of the incoming wave), and exponentially small for the evanescent modes (the higher 
order cavity resonances). Whenever a new cavity mode becomes resonant, the conductance makes a step of universal 
heighth 2e^/h. 

It is well known that multiple scattering by a large density of impurities changes this ballistic transport to a diffusive 
one. In the diffusive regime all intensity is completely randomly distributed over the system. 
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In the intermediate regime of semi-ballistic transport a moderate amount of scatterers is present. On one hand, 
the cavity modes of the pure system are hardly perturbed while on the other hand multiple scattering dominates 
transport in the long direction(s). 

In such systems, interesting effects appear especially near the onset of new resonances. The conductance of a 
quantum wire shows a dip just before a new mode becomes resonant. ||] Certain GaAs/AlGaAs double barrier 

quantum wells show diffusive broadening of their transmission resonance [g| . 

We study the average transmission properties of semi-ballistic systems using the scalar wave approximation in the 
limit of point-scatterers. We assume that neither finite-temperature effects nor Anderson localization are relevant. 

In ref. Q one of the authors discussed a model to explain the transport properties of these systems. In the present 
work we will present the derivations of the results in that article and extend the approach beyond the second order 
Born approximation to the t-matrix. 

The setup of this paper is as follows: In section |l| we will calculate the transmission coefficients of a long, moderately 
disordered waveguide, to supply the derivation of the results presented in [Q. We then apply them to the calculation 
of the conductivity of disordered quantum wires. 

In section [V the transmission of semi-ballistic double barrier structures is examined. Ballistic double-barrier 
structures have a transmission peaked at certain discrete wavevectors. In the semi-ballistic regime there is a broadening 
of these transmission resonances, and transmission to all other resonant channels. 

In section ^ we extend the discussion to include resonance effects of the scatterers that may be induced by the 
geometry. 

In section ^ we discuss the comparison between our model and numerical simulations of an Anderson model in 
the semi-ballistic regime. We close with a summary. 



II. SEMI-BALLISTIC TRANSPORT IN A WAVEGUIDE 



We first consider the case of a moderately disordered rectangular waveguide. To simplify the problem we consider 
scalar waves instead of vector waves, thus neglecting polarization effects. In this way we model the propagation of 
TE modes in electromagnetic waveguides. The same approach also applies to electron propagation in quantum wires 
and, to some extent, sound propagation in long corridors. 

Our geometry has also been chosen as simple as possible. It consists of an infinitely long waveguide with in the middle 
a section of finite length in which a moderate amount of scatterers is present (see fig. |l|). Outside the disordered region 
we have a clean system (perfect leads, for electronic systems) in which ballistic transport occurs. In the disordered 
region semi-ballistic transport occurs. Here we consider (longitudinal) transport inside such a geometry, so that the 
wall potentials can be considered infinitely high. In later sections we consider (tranverse) tunneling through such a 
device and its boundaries. For that apUication the wall potentials must have a finite strength. 



A. The pure system 



Our waveguide consists of four conducting plates, aix — 0,x — dx,y — and y = dx, and it is infinite in the 
z-direction. The conducting plates impose the boundary condition ^' = at their surface, to describe TE waves: 

^iO,y,z)^^P{dx,y,z) = (2.1) 
ip{x, 0, z) — '4'(x, dy, z) — 

For electron propagation this would correspond to infinite wall potentials. We assume that monochromatic waves of 
frequency loq, and free space wavevector fcp propagate through the guide. In the scalar wave approximation this is 
described by the following equation: 

{-V2 - fc2}^(r) ^ 0. (2.2) 

In the absence of impurities the waves travelling through the guide have the form 

V^(a:,y,z) = ^vI;p(p)ei'?P^ (2.3) 
p 

where p — {x,y) is the transversal position and 5'p(p) are the discrete transversal eigenmodes of the system: 
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fplp) = \J^sinip^x)^l^sin{pyy) (2.4) 

Px = —j— Py = —1— (m^j^y > 1 ; m^^y integer) 

Other geometries e.g., cylindrical waveguides, can be described in the same way by substituting the corresponding 
mode wavefunctions for '5. The modes for which < /cq, are 'free', i.e., they can propagate in the z-direction, with 
the mode-dependent wavenumber qp — y/k"^ — p^. The Green's function of this pure system, to be denoted as G*^, is 
the solution of the equation 

{-V2 - fcg}G"(r, r') = ,5(r - r') (2.5) 
with the boundary condidtions (Ejl). To incorporate the boundary conditions in the equation we study the projections 

G° (z, z') = Jdp dp' vI';(p)G°(p; z, pi- z')*p(p') (2.6) 



Since we have translational invariance we can use a Fourier transform to find from ( |2.5| ): 
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The extra term iO ensures convergence of the back-transformation. In real space one has: 



^ ,^ ^ iexpiy^^^^ (2.8) 



B. Adding impurities 

Now suppose a small density of scattering impurities is present at random positions in the region < z < dz- 
Their scattering properties can be expressed in a scattering potential V{r)= — R-i)- We replace the single 

scatterer potential Vg by a (S-potential with scattering strength —u. This is known to be a good approximation for 
electron scattering off a screened charge. The density of scatterers 

n(r) = n{x,y) < z < dz , ^ 

= z<0; z>dz ^ ' 

needs not be homogeneous in the x,y directions, we will assume though that it is independent of z for < z < dz- 



The density of scatterers will be assumed to be so small that the cavity modes (2.4) are still well defined, i.e. the 
scattering mean free path ip must be much larger than the transversal sizes dx,dy. The wave equation for this system 
is 

{-V2 -uYjir- R.)W(r) = P,^j{v) (2.10) 



Here R^ are the random positions of the scatterers, distributed according to a density n(r). 

C. The matrix of a single scatterer 

We first consider the case when only one scatterer is present at position r. The t-niatrix of the point scatterer is 
simply the sum of the Born-series expressing repeated scattering events at the same scatterer: 

<(r) =u + uG{r, r)u + uG{r, r)wG(r, r)u + ... (2-11) 
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Note that since the scatterer is pointhke, there is no momentum dependence and the i-matrix is diagonal in real space. 
Eq. ( 2.12| ) depends on the return Green's function G, which we have not expressed yet. As G is a property of the 



(local) environment of the scatterer it is clear what the physical significance of the ^-matrix is: the t-matrix describes 
the effect of a scatterer in its local environment. In the case where many rando mly positioned scatterers are present, 



the return Green's function G(r,r) in t depends on t itself, which makes eq. (2.12) self-consistent. This means the 
t-matrix is not a property that can be taken from literature. It must be calculated explicitly using the appropiate 
return Green's function in the system under consideration. 



In one dimension {d — 1), equation ( 2.12 ) is well defined. For d > 2 the real part of the return Green's function 
diverges and a reinterpretation is needed. Indeed, for a pure system in three dimensions the divergency is well known 
from the equivalent of Coulomb's law 

G(r, r') = « ^ , „ + ^ (2.13) 

^ ' 47r r - r' Anlr - r'\ At: ^ ' 



As this divergent term appears in the denominator of ( p.l2| ), the t-niatrix vanishes, strictly speaking. This problem 
was discussed for scatterers in free space in ref. We will examine this problem in more detail for constricted 
geometries in section ^ 

For now we will proceed using the simplest way around this problem, known as the second order Born approximation. 
This approximation is commonly used in electronic systems. Indeed, approaches with random potentials that obey 
Gaussian statistics are equivalent to the second order Born approxiamtion. One of our aims is to see whether this is 
still a good aprroximation in cavities. In this approach one has 

tgorn = u iu^Im G(r, r) w iw^Im G(r, r) (2.14) 

The real part u gives rise to a small average potential and will be neglected from here on. This approximation 
maintains the property of scattering but it does not take into account possible resonant behaviour of the scatterer. 
In general, it is a good approximation for very weak scatterers. In section ^ we will turn to the problem of including 
the full t-matrix in our calculations. This will give rise to interesting resonance effects near the subband edge. 

D. The amplitude Green's function and the selfenergy 

The self-energy I](r, r') is defined as the sum of all irreducible scattering events that may be inserted in a Green's 
function line. As the density of impurities is low, we can restrict ourselves to the lowest order approximation to the 
average self-energy, which is diagonal in the space coordinates, S(r,r') = 8(v — r')S(r) with 

I](r) w n(p)tBor„(r) = iu^n{p)\^G{v , r) (2.15) 

The average Green's function G is expressed in the (average) self-energy by the Dyson equation: 

G(r, r') = G°(r, r') -f ^ d3r"G°(r, r")S(r")G(r", r') (2.16) 

To deal with the multiple scattering problem we have to average over impurity positions, weighted with their density 
n(p). Within the second order Born approximation we find an explicit form for the Green's functions Gp of the mode 
P 

In the second order Born approximation I]p=irp, so the resonance width reads 

rp = y"dV?^(p)ImtBor„(p)*p(p) (2.18) 
= y" A^pu^nlyp) lrRG{p,z,p,z)^l{p) 

- J d^puMp)T.^^^P'('^''^)K'(p)K(p) 
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The result for Fp does not depend on z, since after averaging we have translational invariance when we are far away 
from the leads at z = and z = dz- The form (|2.17| ) for G yields in real space 



Gp(., = ■^^MW^j_Pll^Z ^ (2.19) 
2qp + i/ip 



with Qp = Re ^k^ + iVp (2.20) 
and £p = , „^ — (2.21) 

2imyfc^^2qr7jy 

The quantity ip is the mode-dependent elastic mean free path. We find the selfconsistent equation 

Tp = (2.22) 

Tp^E^PP^^P' (2.23) 
p' 

1 -Q 

Vpika) = Rc = = r (2.24) 

2V^-2-p2 + irp |fco'-p' + irpl ^ ' 

Npp, =u^J dV n{p)^l{p)^l, ip) (2.25) 
The number of states (per unit length) in a mode J\fp is 



Afp ^ -Rc Jk^ - p2 + iFp (2.26) 



from which we can see that i^p is proportional to the density of states 

_ TT dA^p 



(2.27) 



E. The Bethe-Salpeter equation 

To describe transport of intensity, electromagnetic energy or the probability of Schrodingers particles through the 
system, we need the averaged intensity Green's function, 

H{r,r') = G(r,r')G*(r,r') (2.28) 

In our model of discrete eigenmodes we consider the projection of the intensity Green's function Hppi. It describes 
the propagation of intensity from mode p to p', and it obeys the following Bethe-Salpeter (BS) equation: 

Hpp,{z, z') = Gp{z, z')G*p{z, z')5pp, + (2.29) 

J dz"Gp(z',z")G;(z',z")x 

p" 

Upp^^Hp.p,{z",z') 

This equation involves the irreducible vertex C/pp'. In our situation it is independent of z for < z < d^, while it 
vanishes in the 'leads' — oo < z < and dz < z < oo. The irreducible vertex can be shown to be the sum of all 
two-particle irreducible diagrams that can be inserted in the intensity Green's function. Two-particle irreducible in 
this context means that the diagrams cannot be split into two separate diagrams by cutting one propagator and one 
complex conjugated propagator line. 

The irreducible vertex is not available in closed form so it must be approximated. We are however not free in 
choosing how to approximate it: the approximation must be consistent with the approximation to S we made earlier. 
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This can be understood as follows: U describes the emission of diffuse intensity by the scatterers, S describes the 
intensity extinction due to the scattering. If the two are not balanced, our description of the system will show gain 
or absorbtion, which is certainly unphysical in the systems we consider here. Flux conservation is guaranteed by the 
Ward-Takahashi identity which can be derived from field theory or by manipulation of diagrams 



i guaranti 

m 0- 



Im j:^=J2Upp'lm Gp-(z,z) 



(2.30) 



If this identity holds, flux is conserved to every order i n the scat terer density while most of our other approximations 
are valid only in leading order in density. Comparing (2.30) to (2.23) shows that 



pp' 



(2.31) 



satisfies the Ward-Takahashi identity. It is, in fact, the ladder vertex, constituting one step in a ladder diagram. It is 
known from transport theory that the ladder diagrams describe diffusive transport. 



F. Solving the transport equation in a waveguide 

Since we cannot solve the Bethe-Salpeter equation analytically we will try to gain as much information as pos- 
sible from approximations. The quantities we are interested in are the average longitudinal intensity transmission 
coefficients Tab, where we use a for the transverse momentum of an incoming cavity mode and b for an outgoing 
mode. 

We consider a wave coming in from z = —oo of the form V'in(r) = a{p) eyip{iqaz) . It will be attenuated upon 
entering the disordered region. This is described by the amplitude Green's function. In the disordered section it gives 
rise to a source intensity 

^p(^) ^W^/^P (2-32) 

We have neglected possible surface reflections as we assume the dispersion relation to be the same inside the disordered 
region as outside. The total intensity present in any mode as a function of position is denoted ^piz). Using ( ^.29 ) it 
can be shown that this quantity obeys the ladder equation 

$p(z) ^Sp+ [ ~ dz'Gp(z, z')G*p{z, z') X (2.33) 
Jo 

E^pp"fp'(^') 
p' 

We can re-express this by expanding the Green's functions 

$p(z) = Sp{z) + (2.34) 

'^dz'e-|^-^'l/^pEc/pp,$p,(^') 
p' 



2rp£p jQ 



where we have inserted the relation 



2vp 



|fco'-p2 + irp| r 



p^p 



Equation ( p.34| ) is a linear system of Fredholm integral equations of the second kind. The solution of this type of 
equation, for oo, is the sum of a homogeneous solution, <I>^, and a special solution which is dependent on 

the source term S^- (In fact, the special solution is defined except for a multiple of the homogeneous solution which 
we can add to it. We will choose such that it remains finite as z ^ oo.) 

By differentiating the system of Fredholm equations we find the following equivalent set of differential equations: 



K(^) = i*p(^) - -FT? E Upp'^'iz) (2.35) 
«p ^p^p p, 

£p$;,(0) = -25p(0)-<i>p(0) (2.36) 

Ep<f'^{d.) = %{d.) (2.37) 
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We will first study the solution to the homogeneou s for m of the system ( ^.35 ) , that is to say, we take S = and 
^ oo. Then we have only the boundary condition ( 2.36| ) at z = . This will then be used to construct to solution 
for finite d,. 



Using eq. ( 2.23 ) the Ward identity ( 2.3C ) can be written as 



(2.38) 



p' 



This implies that the right hand side of eq. (2.35) vanishes if we insert ^p{z) cx i^p. The differential equation therefore 
has a solution of the form 



^^{z) = {zo + z)i^p, forz»£p 



(2.39) 



Near the boundary there will be other terms because of the condition ( ^.36 ). They are related to the non-zero 
eigenvalues of the matrix, so they decay exponentially away from the edges. The asymptotic behaviour of the 
homogeneous solution is characteristic of one-dimensional diffusion: the intensity decreases linearly with z. As 
expressed by the factor i^p in ( ^.39[ ), the intensity is distributed over the modes according to their density of states. 
The shift zq will be calculated further on, when we take in to account the boundaries. 

The special solution for an incoming wave of unit intensity, ?/)'"=e^'?^^5'a in mode a is called In the case of a 
semi-infinite system we can choose it such, that it converges to a constant away from the boundary. The distribution 
over the modes is then given by 



Cg^Vp , Z» £p 



(2.40) 



The coefficient Ca is different for each incoming mode. 

We now examine the behavior of the solutions to (2.35) in terms of the eigenvalues of the matrix of the system. 
There are exponentially growing solutions, exponentially decaying ones and linear -|-constant solutions corresponding 
to the zero eigenvalue of the system. For a semi- infinite system, the exponentially growing solutions will be absent. 
The equation can then be solved formally, yielding 



%{z) = J2 c.^pe"^^' + {a + (3z)Rl 



(2.41) 



where i?p are the right eigenvectors of the system (2.39) and all eigenvalues Xi are positive. The linear plus constant 
term corresponds to the eigenvalue zero of the system, with the right-eigenvector Rp = i^p. 

The boundary condition at z = puts constraints on the coefficients a,/3 and q. For the homogeneous solution 
defined in (2.39) the boundary condition is 



^p<i>;,(o) = $p(o) 



(2.42) 



according to this definition, (3=1 and a = zq which leads to the equation for the coefficients of the homogeneous 
solution c^: 



J2cfrl{Xdp + l) + zoRl = Rlep 



(2.43) 



For the special solution to the problem of a source intensity in channel a the definition (2.40) leads to a = Ca and 
(3 = 0. The resulting equation for the coefficients cf of the special solution reads 



J2 cfi?p(A.^p + 1) + Cai?° = 2(5a,p 



(2.44) 



This equation is very similar to ( 2.43|) , if we take the sum J^a ^a^a--- on both sides of ( ^.44 ) the equations become 
identical and we can conclude 



zo 



(2.45) 



A very useful sum rule can be found from this equation by multiplying (p.44|) by R^ on both sides and summing: 
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^ i?Ocfi?j,(A.£p + 1) + RlC^Rl - 2Rl (2.46) 

i,a 

Taking the inproduct with a vector orthogonal to the leftmost term, but not to R^, we find: 

^ Z^aCa = 2 (2.47) 

a 

To find another useful constant we study the two-particle Green's function H. From the time reversal symmetry of 
the problem we can derive an useful identity: 

Hpp,{z,z') = Hp,piz',z) (2.48) 

If we let z' become large, the Green's function will behave like the homogeneous solution <I>^ near z = b ecau se all 



other contributions are extinguished within a few mean free paths from z' . From the symmetry property (2.48) and 
the behaviour at large z we find the mode distribution at z' must be proportional to Vp: 

lim Hpp^iz, z') = C^vp,<^l{z) (2.49) 

z — >oo 

To find an expression for the coefficient Co we consider a point source in mode a at a large distance zi » ^ from the 
boundary, so that all contributions that correspond to the nonzero eigenvalues of the system will have damped out 
there. It follows we only have to consider the contribution that corresponds to the zero eigenvalue, which yields, by 
considering the jump in the derivative, 

Co^l^Vprpj (2.50) 

Using this coefficient we find an expression for Cpi 

Cp- / dze-^/^PCo Vc/pp'$^,(z) (2.51) 



leading to the interesting relation 



G. Transmission coefficients 

If the sample is finite a certain fraction of the intensity that enters the sample at z = will be transmitted to 
z = dz- We can calculate the transmission coefficients from channel a to channel b for 'optically thick' samples (length 
of many mean free paths) by matching the solution of the ladder equation near both boundaries: 

• For > z but z — dz » £p the solution wil be the sum of the special solution and a multiple of the homogeneous 
solution 

$p(z) = -c$H(z)+ (2.53) 

poo ^ 

J2 dz'i/pp.(z,z')C/p'ae-^/^- 

p, ^0 

• For z < dz but z >> ^p the problem can be considered from z ~ L. There is no incoming intensity and only 
the homogeneous solution will be present, 

<i>p{z) = c^^idz - z) (2.54) 
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In the bulk both solutions have a linear+constant form. This makes it possible to match a (special+homogeneous) 
solution at z = to a (homogeneous) solution at z = dz- 



$^"'''(z) = (a-c(z + zo))i/p (2.55) 

^ ^(^^ _ ^ + ^2.56) 

It follows that 

Ca (2.57) 



(4 + 2zo) 

Since the average Green's function extinguishes in few mean free paths, we do not have to take into account the precise 
behaviour of the intensity at z = for calculating the transmission to z — dz. As usual, since the sample's length is 
many mean free paths, the transmitted fraction of the unscattered intensity is negligible. The intensity transmission 
coefficient for transmission from channel a to channel b is then equal to the intensity ^bidz), 



Tab= / 'dzGb(z,4)Gb(^,4) Vc/bp$p(z) 



{dz + 2zo) 

CaCh 
4{dz + 2zo)ql p 



We have now derived formula (4) of ref . [Q . It holds for the transmission of scalar waves through waveguides and can 
mutatis mutandis be applied to the propagation of EM waves or Schrodinger waves. Below it is used to calculate the 
conductance of a quantum wire. 

III. CONDUCTANCE OF ELECTRONIC SYSTEMS 

We now apply our result to the electronic case of a disordered conducting channel. The Landauer formula gives the 
average zero temperature conductance of a sample of arbitrary dimensions connected to two reservoirs of electrons, 
in terms of the average flux transmission coefficients: 



a.b a,b 



where 2e^//i is the quantum of conduction, the factor 2 comes from spin degeneracy. The flux transmission coefficients 
differ in our case from the intensity transmission coefficients by a factor qb-o/qa-o, where the go's stand for the z- 
wavenumbers of the incoming and outgoing waves outside the disordered region. To a good approximation it holds 
for the propagating modes. 



q.;0^\/ql + i.^^T^ (3.2) 



With these approximations and relation (2.47), we find 
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From this formula, it is easy to see that Ohm's law is valid for the average conductance of samples in the semi- 



baUistic regime. The resistance, defined as the inverse of the average conductance, reads |12 



^ = :r^- + 2^^ (3-4) 



with the conductivity 
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--T^^Erp^^-p (3-5) 

and the 'contact resistances' at z = and z — 

i?c = - (3.6) 
a 

Note, however, that Rc as well as a are complicated functions of fep. 

In ref. an analytic result was obtained for the conductance in a 2D film, with the width rfj, — > oo, in the limit 
of weak disorder: 

where d* — n/kQ is the resonant width, N is the number of open channels and (Tbuik = -ffi^fc/^buik is the Drude 
conductivity of a 3D bulk sample. This result is reproduced here in figure (^, together with the result for a \D 
quantum wire. Both curves have been scaled with the bulk conductivity. 

It is seen that these curves exhibit remarkable drops in the conductivity whenever a new cavity mode becomes 
resonant. These drops are explained mathematically by the density of states that grows large near the subband 
bottom, thus causing the second order Born i-matrix to become larger. This expresses more efficient scattering and, 
therefore, less conductance. The physical explanation is: When the 'new' mode is not yet resonant it does not yet 
contribute to conduction. There is scattering to this mode however, and the scattered waves interfere destructively 
with the waves present in the other modes. In section ^ we will study the analogue of this effect for the full i-matrix 
of the pointscatterers. 

A different approach for calculating the conductivity was followed by Surke and Wilke in ref These authors 
calculate the conductivity directly from the Kubo formula and derive expressions different from ours. Their results 
involve the average scattering time, rather than the inverse of the average scattering rate (Fp in our work). It therefore 
seems to us, that the latter results are unphysical. Indeed, one can consider cases where the scatterer density goes 
to zero locally, such that the average scattering time diverges. Then the resulting prediction for the conductivity 
diverges, while our result (|3.5|) is quite insensitive to such limits, as it should be. 



IV. TRANSPORT THROUGH A DOUBLE BARRIER STRUCTURE 

In section transmission coefficients were derived for transport of wave intensity along the length of a waveguide. 
In some systems, like the Fabry-Perot Interferometer (FPI), or its electronical analogue, the double barrier quantum 
well (DBQW), transport occurs in the transversal direction due to 'tunneling' through the barriers. 

In the absence of random scattering these devices transmit only waves for which the perpendicular component of 
the wavevector is resonant with the cavity. The linewidth is very small, usually the Q-factor of these devices is several 
thousands. The pure FPI transmits a light beam that meets the resonance condition without changing its direction. 

The devices we are interested in contain a small density of impurities, such that the width d^ of the device is still 
much smaller than the mean free path, but the width multiplied by the Q factor is much larger than the mean free 
path. The problem of multiple scattering in such devices was first considered on a fundamental level by Berkovits 
and Feng [[l4[ . Their approach is valid for the situattion of one resonant mode, well away from the onset of further 
cavity resonances. 

The behaviour near resonances in the multimode situation was later discussed by one of us jj]. In the present 
section we will discuss the derivation of these results. It will be seen that multiple elastic scattering will broaden 
the resonance linewidth and cause transmission of the energy into all available outgoing channels, independent of the 
incoming channel. 



A. Double barrier system in one dimension 

First, for simplicity, we will consider a one-dimensional double barrier quantum well, like in fig. (^, with the 
potential between the barriers equal to that outside. We will describe these barriers, which are imperfect mirrors, by 
strong (5-function potentials. Their strengths are allowed to have an imaginary part to model absorbtion of waves by 
the mirrors. 
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V{x) = Vj^5{x) + VnS{x - d^) (4.1) 

The; intensity transmission, reflection and absorbtion coefficients of the single barriers follow from requiring conti- 
nuity of the wavefunction across the barrier: 



^ |y + 2ifcoP ^ |l/ + 2ifco|2 ^^'^^ 

If V is purely real, there is no absorbtion and the identity T + i? = 1 holds. If Im F > the model describes a barrier 
with nonzero absorbtion coefficient A: 

„ _ 4fcoIm V 

A description with absorption is realistic in the optical case, where the absorbtion coefficient of a mirror may be 
larger than the transmission coefficient. In the electronic case there should be no absorbtion, but there could be loss 
of phase coherence which can lead to similar effects. 

Now we look at the Green's function inside the structure, G{x,x') for < x,x' < dx- By requiring continuity at 
both barriers we can see it is, neglecting terms of order l/V"^, 



KO Sm^fcod^J + VL-2iko + VH-2iko 



where J is a (complex) phase shift which becomes small near resonances. In the ideal system {V real and very large) 
these resonances occur when ko = rnr/dx- If the linewidth of the structure is much smaller than the mode separation, 
which is the case for high-Q cavities, we can make a Lorentzian approximation to the lineshape and express the 
propagator as a sum over mode wavefunctions: 



G{x,x')= '^pA^)GpAko)^pA^') (4-5) 

„ _ nTT 
~ da, 

involving the mode wavefunctions 

(a;) = \ -T- sm{pxX + Sp^ ) (4.6) 

V 

and the mode-propagator, (upon neglecting small shifts of the resonance positions) 

Gp^ (ko) = J^ypnre (4-7) 

The linewidth TP""^® is due to absorbtion and transmission of the barriers: 

^pT = ^{AL+TL + An + Tr) (4.8) 

where L, R stands for left and right, respectively. This expression of the linewidth in terms of transmission and 
absorbtion coefficients is valid for general barrier potentials. 

A beam coming in from x = — oo of the form exp(iA;oa;), with ko close to some resonant Px will give rise to a 
standing wave inside the structure. The amplitude of this wave is found by matching at the boundaries, which gives, 
approximating for ko near Px : 

To describe realistic situations we need to expand the simple ID model to a three dimensional model with barriers 
('mirrors') a.t x = 0,x = dx- There are three regimes for this generalization: 

• A structure which is confined in both y and z dimensions, like a double barrier structure inside a waveguide. 

• A structure which is confined in the y direction but very large in z dimension, for example a double barrier 
structure inside a 2D electron gas. 
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• A structure which is very large in both y and z dimensions, hke an optical FPI or a double barrier quantum 
well inside a three dimensional semiconductor sample. 

The first case is the simplest to describe, we will analyze it first and then take the continuum limit in two directions. 
Results for the pure case follow from simply replacing the ID mode eigenfunctions with the 3D eigenfunctions and 
replacing the mode energy by = p'^ + + Px 

G(r,r')= *p(r)Gp(fco)*p(r') (4.10) 



*p(i") = y ^ ^ ^ sm{xpx) sm{ypy) sin(zp^) (4.11) 

We can see from the diagonality of the Green function that only waves which have a resonant with the structure will 
be transmitted into the same y, z mode of the waveguide surrounding the structure, i.e. the interferometer selectively 
transmits waves, it does not change their direction. 

In the case where scatterers are present, there will be a self-energy TP due to disorder in addition to irP"'°. One 
has 

Ep = SP">-^ + w iPP"'-" + (4.13) 

We are interested in the regime where the latter has a much greater imaginary part, r^^rP"''°. This is the case for a 
cavity with a very high intrinsic Q-factor and a low density of scatterers. Let us consider the transmission of a plane 
wave coming in from x = —oo^ with wavevector a. (Note that in this section we consider transmission through the 
cavity boundaries whereas in the previous section transmission was within them, therefore the incoming nbeam is now 
in the x-direction.) The momentum in the y and z directions is quantized while the momentum in the x direction 



outside the structure is determined by the total wavenumber. 



= e^axx^_±^ sin{ayy) sin(a^z) (4.14) 
This incoming wave gives rise to an intensity inside the structure of 

- E f^J^p |G'pn^p(r)P<5p,.a,<5p.,a. (4.15) 
In the case no scatterers are present, the intensity of the transmitted wave will be 

= E §T^T^\G^\H,^^aJ,^,a. (4.16) 

Px ^ 

At resonance this transmission is unity for the case of symmetric, non-absorbing barriers. In practice, only one of the 
modes will contribute to this transmission. 

In the case when scatterers are present, intensity will build up in the cavity in essentially the same way, due to the 



condition I » d^, but the dominant decay of the intensity will be due to scattering to other modes. Hence, (4.15) 
is still valid in this case, but we must take the Green function with the disorder term Tp in the self-energy. The 
incoming intensity will now be the source term for a BS equation which reads 



*P = IGpP <^ -^T^Sp^,a,6p^,a. + E \ (4.17) 




$p is the intensity in mode p, t/pp' is the irreducible vertex. Since there are only discrete modes now, this is just a 
matrix inversion problem: 
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Im G 



^t/pp' ) <i>p' = |Gp 



2 rpL 



Tp ^Py,aySp^^a^ 



(4.18) 



The matrix in the brackets has one eigenvalue close to zero which will dominate the solution. To see this we multiply 
by Fp and sum over p, using (2.30) and (4.13) to find 



(4.19) 



We expect the intensity in each mode to be proportional to the density of states (this approximation is very good 
if the scattering rate is much greater than the cavity loss rate). Also, if the in comin g wave is close to one resonant 
mode we can neglect contributions of the other modes. This allows us to solve ( 4.18 ) and find 



*n = 



TplGpp 



' Ep'rp"'=rp,|Gp,p 

The transmission from channel a to b is then given by 



(4.20) 



(4.21) 



In case the double barrier structure is large {dx » £) and uniform in one or both of the {x, y) dimensions, the sum 
in the denominator can be simplified since the resonance width ^p^,py,p^ will no longer depend on the corresponding 
wavenumber. For the case the structure is large in the z dimension we find 



ab 



a^ri^ralGapT^^rblGb 



dld.Y. 



(4.22) 



with the definition of Vp given in ( 2.24 ). This is eq. (7) of ref. In a more restricted version with one resonant 
mode it was already derived in ref. |14| ; note, however, that in that work the prefactor misses a factor 2. 

In the case of a double barrier structure that is much larger than the bulk mean free path in two dimensions. 



[dx, dy >> £), the denominator of (4.22) can be simplified to 



Mldyd, 



Eppurc 
Pre 

Pi 



(2- 



■ arctan ■ 



(4.23) 



/^From ( 4.22 ), we see that scattering occurs into all free cavity modes. For a FPI {dy oo) with incoming plane 
wave, scattering inside the barrier region causes an equal distribution over all directions. Therefore, the transmitted 
intensity does not depend on the angle of incidence in the (y, z) plane. It does depend on the angle with the x axis 
as the mode quantization implies that only certain angles are transmitted. This causes the well known appearance of 
fringes in the transmission pattern. 

When a new cavity mode is just resonant, the return Green's function in the cavity becomes large, as is discussed 
in section VF. This quenches scattering in the other subbands, so the device will appear less disordered. This is one 



of the reasons why a Fabry-Perot etalon gives best results when used near perpendicular incidence. We should point 
out, however, that this effect is less pronounced for a two-dimensional structure such as a Fabry-Perot interferometer 
than for a one-dimensional structure such as a wave guide, as the divergence of the return Green's function is only 
logarithmic. 



B. Electronic conductivity through a disordered double barrier structure 

In results were also derived for electronic double barrier quantum wells. The conductance of such a device is 
again given by the Landauer formula, 

ab 
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In an experiment conducted by Gueret et al, the barriers consist of two separated layers of aluminium gallium 
arsenide (AlGaAs) in a gallium arsenide (GaAs) sample. 

The thickness of the barriers was varied from 7.5 nm to 31 nm. Hence, the transmission coefficients, which depend 
exponentially on the barrier thickness, varied over nine orders of magnitude, while the width of the resonance peak 
remained a few mV. Elastic scattering from roughness at the interface between the GaAs and the AlGaAs is believed 
to cause this resonance broadening. 

For a 'pure' double barrier structure, the linewidth is expected to scale linearly with the transmission coefficients: 

r = T'^^^iAL+TL + AR + Tr) (4.25) 

The total transmission for a resonant cavity is linear in the transmission coefficients, as it is proportional to TlTrF^. 
In the case of a small amount of disorder, for a cavity at resonance, the transmission coefficients can be seen from 



(4.22) to be proportional to TlTrF^, as in the pure case. The linewidth, however, is determined by F^, which can 
be much larger. It was calculated in |j] that this multiple scattering effect can explain the order of magnitude of the 
observed resonance broadening. 



V. STRONG SCATTERING: THE FULL BORN SERIES 



The results of section ^ have been derived using the second order Born approximation ( p.l4| ) to the t-matrix. 
This approximation essentially describes weak scatterers that have no internal resonances near the wavelength of the 
incident waves. Realistic systems, however, are not extremely long. If a reasonable amount of scattering is to occur, 
one must have moderately strong scatterers. But for such scatterers the multiple scatterings from the walls will not 
at all be negligible. Rather than approximating the Born series by the first two terms one has to sum the full series. 
For point scatterers in continuum space this requires introducing either an extension of wavefunction space [ p^ or 
a regularization of the return Greens function. A very simple solution to this problem was given in ref. In the 
present section we shall apply those ideas to cavity systems. 

Scattering of Schrodinger waves in narrow systems has been studied by Bagwell ||l^. He considers both finite size 
scatterers and point scatterers. In the latter case, the infinities that occur due to the small distance behavior are 
avoided by considering only a finite number of evanescent modes. As a result, the predictions of the model then 
depend on the number of modes considered. Chu and Sorbello [Q investigated scattering properties of finite size 
scatterers inside a cavity by means of the imaging method. We shall follow a different approach that is more suitable 
to the limit of point scatterers, and then proceed to multiple scattering situations. 



A. Point scatterer near a mirror 



Boundary conditions can have a strong effect on the properties of scattering, even for point scatterers. This can 
already be seen from studying the following simple case: Consider a point scatterer for Schrodinger waves with 'bare' 
scattering lenghth m at a (small but finite) distance dw from a perfect plane wall. This imposes the mirror boundary 
condition ^' = on the surface z — Q. The eigenfunctions of the pure system (without the scatterer) are 

*(p.,Py,P.) = V2ie'''^^e'P«^sin(p,0) (5.1) 

with 

Vl+Pl+Pl^kl (5.2) 

The return Green function is found by Fourier transformation: 

G(ro, ro) = G°(ro, ro) - G"(ro, r^) (5.3) 

Where = {r1,ry,—r1) is the mirror image of Tq, and G° is the free space propagator. We distinguish a 'direct' and 
a 'mirror' term. When the distance between the scatterer and the mirror is much larger than the effective size of the 
scatterer, (d^, 3> A^^), the mirror term has a regular behaviour. 



We then need to regularize only the 'direct' term in (5.3), the regularized 'mirror' term will not be much different 



from the unregularized form. This enables us to apply the regularization scheme in |9| to the direct term only, to find 
corrections due to the scattering between the impurity and its mirror image. It gives: 
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(5.4) 



However, when the distance between the mirror and the impurity is of the order of the effective scatterer size, this 
approach becomes ill behaved and we need to regularize both terms in the return Green's function. We do this by 
applying a regulator of the form A^/(A^ + P^) to the Fourier integrals that are used to define the real space form 
of the Green's function. This regulator is well known in quantum field theory. It cuts off momenta beyond p = A, 
showing that 1/A is the effective size of the scatterer. In the free space situation we thus get 



1 



A2 



(27r)3 p2 - 
^i\r-v'\ko 

47r|r - r'l 



iO A2 



r'|A 



,ip • (r - r') 



A2 



47r r - 



A2 



Now we can take the limit r' 



r. As expected, the result is finite, and reads 
ifco + A A2 



G"(r, r) 



47r A2 



ftp 



(5.5) 



(5.6) 



The factor A2/(A2 + fc2) jg near unity for scatterers which are much smaller than the wavelength. We will omit it 
from now on. Then (5_^) coincides with the result of ref 1^: 



G"(r,r) 



ifco + A 
47r 



The return Green's function (5.4) now reads 



ifco + A e2*'=^t"^- - e 



47r 



8ndy 



(5.7) 



(5.8) 



and reduces to eq. (5.4) when exp(— 2Ad„) <C 1. Now we use this regularized form of the return Green's function to 
compute the t-matrix 



A + ifco 



g2idu,fco _ g-2d^A' 
87rd„, 



(5.9) 



Following we state that by definition, a position dependent resonance occurs when the t-matrix is purely imaginary. 
(This need not always coincide with a sharp peak in the absolute value of the t matrix) . The positions of the resonances 
are located at fco -values where 



cos(2(itufco) ~ Sird^tU + 2d„,A + exp(— 2(i^A) 



(5.10) 



In case of Schrodinger waves, u is the potential strength and does not depend on fco. This implies that for large dw 
there is no resonance. Indeed, in free space {d^ = oo) it is known that point scatterers for electrons do not have a 
resonance. However, eq. ( |3.10 ) shows that an infinity of resonances occur when the scatterer is attractive {u > 0) 
and is located close enough to the wall. 

For scatterering of scalar classical waves, the analogue of the above position dependent resonance can be found by 
replacing the scattering strength u by AnakQ, where a is the polarizability of the point scatterer The expression 
for the i-matrix becomes 



t = 



1 



'a + ifco e^i^'^fco _ e-2rf«)A^ 



Sird^i 



It has resonances at 



cos(2dt„fco) = 2dw [ - A ) + e 



-2d,„A 



For du, — > cx) this yields the resonance condition kg — = 1/ Aa of ref. 
may occur. 



(5.11) 



(5.12) 

For large but finite dw many resonances 
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B. Pointscatterer in a Waveguide 



N ext, consider a scatterer in an infinitely long channel. The Green's function for the clean channel is calculated in 
(pTq). It is expressed as a sum of mode propagators 



G(r,r') = V*p(p)Gp(z~z')*p(p') 



p 



(5.13) 



The mode quantum numbers Px,Py are the discrete transversal wavenumbers. 

For systems on a lattice the sum over p converges as the distance between adjacent lattice points is a natural cutoff. 

Even in a confined geometr y of c ontinuum space we have the problem that the real part of the return Green's 
function diverges. In equation (5.13) this occurs since for large |p| the terms in the infinite series tend to zero as 1/p, 
leading again to a linear divergency of the sum. In order to define the i-matrix of the point scatterer in a meaningful 
way, we thus still need to regularize the return Green's function. 

The approach followed by Bagwell is to take into account only a finite number of evanescent modes in equation 
(5.13). As the cutoff has not been related to physical properties of the scatterer, this approach is a bit unsatisfactory 
and unnatural. 

We can regularize the sum ( ^.13 ) term by term in a way similar to (5.5), using the cutoff function / [K^ + 'p'^ + q^), 
where p is now the discrete transversal wavenumber and where q is the continuous wavenumber in the z-direction. 



Gp(z,z') 



dq 



1 



27r g2 



q(z-z') 



-p2|z- 



2\/A^ 



A2 



A2 



(5.14) 
(5.15) 



We will consider the real space case again, omitting the factor A? /{h?' 
function will then be given by 



Sp) from now on. The return Green's 



G(r,r)^^*2(p) 



1 



/r2 



2yA^ 



(5.16) 



For large |p| the summand is of order p so the series (|5.13| ) converges 



For the special case where the energy approaches one of the subband bottoms p^ the series contains one divergent 
term. As we shall see in detail below, this divergency will be regularized when the system is finite and there is 
conductive behavior near the boundaries. 

The t- matrix of the point scatterer (in a channel without further randomness) now reads 



(p) 



2vV 



1 



fr2 



1 



2V^ 



(5.17) 



If the scatterer strength is positive, that is to say, for an attractive scatterer potential, the t matrix can have a pole. It 
occurs before the onset of conduction modes, thus for fco < I Pol: where po is the label of the first mode. This bound 
state will have an energy between < fcp < Pp, the energy of the first subband bottom, provided that 



< u < 



1 



1 



n -1 



L P 



2v/p' - Sp 2v/A^T^ 



(5.18) 



Thus, the bound state already occurs for rather weak scatterers, which is an indication of the failure of low-order 
approximations, such as the second order Born approximation. On the other hand, a strong attractive scatterer will 
have its bound state at negative energy, i? = fcg < 0. In an optical situation such a state has no meaning. In an 
electronic system it may occur, however. If the Fermi energy is above the bound state energy, the state will always 
be occupied and it will not contribute to conductance. Due to the divergency of the return Green's function at the 
bottom of the first subband, the bound state will always occur below the first subband bottom. The weaker the 
scatterer is, the closer the bound state energy will move towards the first subband bottom. Also note the explicit 
depencence on the transversal (x, y) position of the scatterer. For very weak scatterers, however, the resonance will 
be indistinguishable from the subband resonance and the second order Born approximation will yield good results. 
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C. Resonant Tunneling 



Unexpected peaks have been observed in the transmission of quantum dots, narrow constrictions of a 2D electron 
gas . Before any of the channels of the dots were open, a transmission peaks of order unity were observed. They 
were attributed to tunneling through a bound state of a single impurity located (by accident) in the channel. We 
model such a narrow junction between two broad reservoirs where conduction takes place by an infinite channel. We 
assume that the electrostatic potential equals Vb for < z < L and vanishes elsewhere. This implies that for z < 
and for z > L conduction occurs, while the 'junction' < z < L acts as a barrier as long as /cq < Vg. 

In the true ID case this problem can easily be solved exactly, yielding the amplitude transmission coefficient 

8ifco/^"e-^ 

^ (2K-u)(fco + i«;)2-ttyo[e-2'=^«+c-2«(i-^o)] + (2K + tt)(A:o-iK)2c-2'=-f' ^ ' ' 

The intensity transmission coefficient is the absolute square of the amplitude transmission coefficient C. It reads, 
up to an exponentially small term in the denominator, 

T= 64fcg.^e-^-^ 

4fc2«;2(2K_u)2 + [(2^_^)(^2_^2)_.yVb(e-2«Zo+e-2«(i-^o))]2 ' 

Since the potential is equal on both sides of the barrier, and the incoming and outgoing modes are the same the 
flux transmission coefficient T equals the intensity transmission coefficient T. When = = Vq — u'^/4:, there is a 
resonance peak in the transmission. Since then k — u/2, its height is 

T„.ax = r(fc.) = v;^{l + coshu{L-2Z,)) ^^-^^^ 

^From this equation we can see that the tunneling transmission is maximal if the scatterer is positioned in the centre 
of the barrier, Zq = L/2. Then, the transmission coefficient is no longer exponentially small. Independent of the 
length of the barrier, its maximal value is 

{Vq - uy A)u^ _ ^ (V^-u^/2 
V6 



This phenomenon is called resonant tunneling. We see that T is unity if — 2Vo. This occurs when the resonance 
energy is exactly half the barrier energy, k^ — Vo/2. We will now show how this feature of a truly ID system persists 
in the quasi- ID waveguide. 



D. Resonant tunneling in a waveguide 

The situation for impurity tunneling through a barrier in a waveguide is analogous to the ID case if we can assume 
that the lowest subband is nondegenerate, i.e. there is a lowest subband pq and all other subbands are located at a 
higher energy. In that case, we can ignore tunneling through the higher subbands. We do need to include all subbands 



in the calculation of the return Green's function, as well as reflection terms for the lowest subband. ^From (5.15), the 
return Green's function becomes 

G(po, Zo, po, Zo) - J2 *p(P)Gp(0) (5-23) 
p 

+Gp„(-Zo)rpGp„(^o) + Gp„{L - Zo)rpGp„{-L + Zo) 

We will use the abbreviations Kp for the imaginary wavenumber inside the structure, qp for the wavenumber outside 
the structure, rp for the internal (amplitude) reflection coefficient at the edge of the barrier. 



= . /p2 + _ ^2 ^ ^ (5 24) 



gp = V ^0 ~ 9 = Ipo (5-25) 

gp - i^p 
' gp + iKp 



= -2^p '^P (5.26) 
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Thus, we can calculate the return Green's function, including the reflection terms 

(5.28) 

Neglecting the direct tunnelhng term we can calculate the transmission coefficient for a barrier with a scatterer at 
position (po, Zq). 

T = . . lefjfi^ 

{2-uR)^K^Vo^ - Ku{2-uR)VQ{q^-K^)[c~^''^''+e-^''(^-Z"')]+uW^[c-^^Zo+c-^'^(L-z„)^2 



with the abbreviations 



R^Y. KiPo) i , , , \, - , , I (5.30) 

w = ^.vi/2^ (5.31) 



The resonance condition now is uR — 2, where the peak flux transmission coefficient becomes 



"^'^ V^(l + cosh2K,(L - 2Zo)) 
Again the maximum peak height occurs when the scatterer is located at the middle of the barrier 

%na^[L/2\ = M^K^/Vq = 1 - 



(5.32) 



Vo 

It depends on the transversal scatterer position and on the scatterer strength only through the shift of the resonance 
position. Its maximal value is again unity. 

To calculate the width of the resonance peak of ( ^.29| ). To do this we will need some derivatives: 

d kf) 

— Kp = (5.33) 

oka Kp 

i-^R-koY.HiPo)^^^HM^^ (5.35) 

The approximation for the derivative of R is good if Kp„ is much smaller than the other decay constants, which it is 
in a usual tunneling situation. 



The denominator of ( 5.29 ) consists of three terms, the first two are zero at resonance. Half height occurs approxi- 
mately when the sum of those terms equals the third term. The first term of the denominator is of second order in 
Afc = A:o — /c*, the second term is of first order, but has a much smaller coefficient, so we will have to include both 
terms to determine the peak width. We solve 

2 2a^ /s: 
to find the result 



-2kZ, 



(Afc)FWHM = vV+^ (5.37) 



-2kZ, 



'0 _j_(-,-2K(L-Zo) 



w ; (5.38) 

ko 

So the peak full width at half maximum, (Afc)FWHM, is exponentially small even if the scatterer is positioned near 
the center. Only through the resonance energy, determined by uR — 2, it depends on the {x, y) coordinates of the 
scatterer or its scattering strength u. This means that of the three quantities that are easily accessible by experiment, 
the resonance position, heighth and width, only the resonance position contains information about the scattering 
strength and (a;, y) position. 
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E. Multiple scattering 



Suppose that in the region < z < there is an a priori nonuniform density n{x,y) of randomly placed, possibly 
resonant impurities. To first approximation, this causes the average subband Green's functions to obtain a self-energy 
term 



Ep= / d'p\^pip)\'n{p)t{p) (5.39) 



The self-consistent t matrix can now be calculated by inserting the self-energy E in ( 5.17 ). In order to describe 
diffusion , we need to define the ladder vertex U, which must obey a local conservation law (Ward identity) of the 
form of ( p!30| ) 



ImSp = Upp'lmGp, = ^ Upp^i^p (5.40) 
p' p' 



This identity ensures flux conservation. We define the ladder vertex U as 

C/pp, EE I d'p^l{p)^l,ip)n{p)tip)t*ip) (5.41) 
Then the Ward identity ( 5.40|) indeed is satisfied, since 



Where we have written G for G(r,r). With the vertex function defined in ( ^.41 ) we can reproduce all results from 
section |l| if we make the following modifications: 

• Wherever the expression /cq — + iVp is used, we need to replace it by /cq — + Sp, thereby incorporating 
the real part of the self-energy. 

• When the imaginary part Fp is used outside the propagator denominator, we should replace it by Im Ep. 

• The matrix A'pp' , the second order Born approximation to the irreducible vertex, should be replaced by the full 



i-matrix approximation C/pp' defined in (5.41) 



To find the correct self-energy we need to solve t{p) and Gp simultaneously, in general we will need to do this 
numerical ly. A s the equations are well-behaved in a physical situation, an iterative improvement method can be used. 



In section VI we will discuss an approach for a lattice model where a direct simulation of the scattering problem was 
performed. 

The drops in the conductivity that can again be seen in the case of attractive scatterer, shortly before the opening 
of a new subband can be explained qualitatively as single scatterer resonances. The real part of the Green's function 
of the not-yet opened subband starts to diverge, which means it has to cross a resonant value, rendering the t- 
matrices of the individual scatterers imaginary and very large. Thus, the scatterers scatter more eficiently at these 
energies, thereby decreasing the conductance. As the resonance is position-dependent, there will not be a sharp drop 
in conductivity, but instead a smooth decrease. When the new subband opens, the density of states effects will cause 
a sharp drop in conductivity, which is already present in the second order Born approximation. 

For repulsive scatterers no resonance exists and the conductivity will rise with increasing energy until the opening 
of the nest subband, at which time it will drop again. 



F. Subband bottom transparency 

For the case of a single scatterer in a channel an interesting effect was observed by Chu and Sorbello and by 
Bagwell ||l^: when the energy is close to the bottom of a subband, the return Green's function is very large so that 
the t-matrix is small. As a result, there is almost no scattering. 

For the multiple scattering situation the question that arises immediately is: will the system become optically thin, 
so that ballistic transport determines the transmission coefficients? 

Let the energy approach a subband bottom feg « p,^. The return Green's function Gp, of thi s mo de will become 



large due to its quasi 1 — D square root divergency. The corresponding term in the denominator of (5.17) will dominate 
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the other terms, except at the nodes of the corresponding subband wave function. If this wavefunction has no nodes, 
or the scatterer density is zero near those nodes, as can be the case for a waveguide with disorder only near the edges, 
the self energy of the corresponding mode can be solved approximately: 



dVl*P.(p)l'ri(p). ^ 



*P.(p)PGp.(0) 



= 2in,„^A:2-p2 + Ep, (5.43) 
where Uid is the ID scatterer density 

n,n = J n{p) (5.44) 



We see that (5.43) does not contain the scatterer strength u. It is valid for large values of u, as for v ery s mall scattering 



strength the term u ^ will not be neglegible compared to the return Green's function. Squaring (5.43) and requiring 
that the imaginary part of the self-energy stays positive, we find 



Ep. = -2n2^ + 2ini^^fc2 ^2 (5.45) 

This expression is of second order in n{p). Although our theory is a first order approximation in nt, we can still use 
it because the t matrix is so small near the subband bottom. This is the multiple scattering analogue of the subband 
bottom transparency effect noted by Chu and Sorbello |^] and by Bagwell for the case of a single scatterer: when 
the energy is at the subband bottom there is very little scattering, and the system may even become optically thin, 
in which case ballistic transport determines the transmission coefhcients. We can calculate the behaviour of the mean 
free path, which is, unlike the selfenergy, usually easily experimentally accessible. 



}-i = 2Im Jfc2 _ p2 _ 2 n?^ + 2ini„^/c2-p2-n2^ (5.46) 



= 2Im (^i n,^ + ^Jki-pl- j (5.47) 

This gives immediately 

^ (5.48) 



If there are lower lying, already propagating modes, their mean free path can be much longer than this. Let's first 
consider the case where the scatterer density n(x,y) is zero near the nodes of the mode wavefunction 'I'p.,. Then we 
can find the self-energy of the other modes by 

If we estimate the value of the integral to be roughly unity, we find for the mean free paths of the other modes 



« 2Im y fc2 _ p2 _ 2n2^ + 2in,„^kl - p2 - n?^ (5.50) 
Since — p'^ is much larger than the other terms we can neglect terms of order (/cq — P*)/(fco ~ P^) to find 



fp}-i«2n,.^^M^^L^ (5.51) 
V fco - P^ 
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Where « fcg is the newly opened mode and << f^'' ^^'^ already open modes. It is seen that this expression 
vanishes at the exact subband bottom, so the mean free path diverges and the transmission becomes ballistic in every 
mode except the newly opened one. 

In practice, the scatterer density will not in general be exactly zero near the nodes of the wavefunctions. The 
t -mat rix may be appreciable near the zeroes of the wavefunction 5'p_^ of the 'new' mode, so that the approximation 



( p. 50]) cannot be made. This will have non-neglegible effects on the behaviour of the other modes that do not have 
the same zeroes. The corrections to the self-energy of the propagating modes that arise will be of order 1/-^ Gp^ , 
which is very small near the subband bottom so the resulting mean free path may, with these corrections, still become 
much larger than in the middle of the band. If it approaches the order of the system length, we can no longer use the 
results derived for optically thick systems. Thus, a sample may cross over from the semi-ballistic regime of multiple 
scattering to the ballistic regime of 'low-order' scattering when the energy is near a subband bottom. 

VI. NUMERICAL SIMULATION 

In order to look numerically at the effect of scatterers in a wave guide we use a model which contains the essential 
physics of the continuum Schrodinger equation whilst being ameanable to this kind of analysis. One such model is 
the Anderson Hamiltonian on a square lattice. 

{-\nx,ny > e„^^ny < n-x^riyl 

,ny 

+ Irix + l,ny >< rix^nyl + \nx - l,ny >< nx,ny\ + \nx,ny + 1 X nx,ny\ -\- \nx,ny - 1 >< nx,ny\} 

By setting the hopping elements to unity we define the unit of energy. The model has an upper and lower limit of 
energy E = [—4, 4] for the ordered case en^,n„ = 0. The electron like band structure we are interested in starts at the 
band edge closest to E = 4. The sign of the disorder parameter tn^,ny has been chosen such that for positive values the 
scatterers are attractive and negative values repulsive. The conductance of this system can be calculated by simply 
transforming the Hamiltonian above into a set of transfer matrices, one for each layer, and multiplying them together 
to find the transfer matrix for the whole system. However, in the strong scattering regime where the conductance is 
very small this method is unstable since half the eigenvalues grow exponentially and half decay. We therefore prefer 
to calculate the conductance from the Greens function which either decays or grows. The method we use to calculate 
the Greens function was introduced by Lee and Fisher and given for the general case by MacKinnon It 
consists of using Dysons equation to calculate the required elements iterativley by adding sucessive disordered layers. 
The general form for the relation between the Greens function elements and the transmission coefficients has been 
given by Baranger and Stone [ pO[ and we use the Landauer-Buttiker formalism to calculate the conductance from 
these coefficients. 

One of the major differences between the continuum model and the Anderson model is in the number of evanescent 
modes in a complete basis at the Fermi energy. The continuum model has an infinite number whereas the Anderson 
model only has as many modes as there are lattice sites across the wire. We can easily accomodate our continuum 
theory to this band structure by taking an appropriate momentum cutoff function. In the simulations we present in 
this paper we have been interested in disordered systems which are long in comparison to the mean free path but 
much narrower. In this limit only the evanescent modes which have any appreciable probability amplitude between 
scattering events need be included. This makes the Anderson model an ideal candidate for this calculation. 

Comparison between these numerical simulations and our continuum model can be made without any free param- 
eters since the scatterer size is fixed by the requirement that each scatterer should occupy the volume of one lattice 
cell. 

In figure ^ we compare the semi-ballistic theory (full line) with one realization of the Anderson model (dotted 
line). The sample configuration is that of a long, weakly disordered wire. In figure ||the same is shown for repulsive 
scatterers. Both graphs also show the same curve calculated in the second order Born approximation, in which there 
is no difference between attractive and repulsive scatterers. The conductance is shown as a function of carrier energy, 
where the units of energy are related to the band structure of the lattice (—4 = zero energy, band center). It is 
seen that the subband maxima have a round shape for attractive scatterers while they are quite sharp for repulsive 
scatterers. Also, the conductance is lower for attractive scatterers. This can be explained as an effect of the scattering 
resonances that exist for attractive scatterers and are described in section ^ 

The shape of the conductance graph is seen to be dependent on the sign of the scattering strength, its magnitude 
and the scatterer density. Depending on these parameters the conductance trace can be either smooth, sawtooth-like 
(fig or have steps (fig |l^. These steps can be (much) smaller than the universal height 2e^//i of steps in the 
conductance of ballistic devices, although they occur at the same energies. 
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It is also seen in these graphs that the second order Born approximation gives a conductance that is much lower 
than the full t-matrix, although it takes into account less scattering diagrams. This effect is shown for two different 
scattering strengths in fig. The explanation is that contributions of higher order scattering processes in the full 
t-matrix have the opposite sign of the lowest order contribution. 

In figures ^ and |9| conductance of shorter samples with a relatively high density of scatterers are shown. The good 
agreement between our theory and the numerical date (average over 10 samples) demonstrates that our model, which 
assumes pointlike scatterers is still good at finite volume fractions. 

In figure ^ we see the contact resistance term i?c and the extrapolation length zq, which gives rise to this contact 
resistance. It is clearly seen that zq is highly nonuniform and varies within each subband with a small dip attributed 
to scattering resonances (the return propagator being small and real there), a peak owing to a weak form of subband 
bottom transparency (the propagator being large and real, so that scatterers cannot be resonant and scattering is 
slightly suppressed) and a steep drop caused by the sudden increase in density of states as the propagator crosses over 
from large and real to large and imaginary. Note that the contact resistance Rc has a much smoother behaviour. It 
occurs since the conductivity a essentially varies in the same way as zq 

The extrapolation length zq is a weighted average of the mode mean free paths ( 2.45 ) so it is of the order of one 
mean free path. Since our theory presently does not include quantum localization effects, it does not remain valid for 
samples which are many mean free paths in length. Therefore, in the regime where the semi-ballistic theory is valid, 
the influence of the contact resistance term in eq. (3.4) will always be important. 

Figure |^ shows an example of a wire which is so long that localization effects can no longer be neglected 

Nikolic and MacKinnon pl[ discuss the conduction of quantum wires with a large amount of boundary roughness 
and a number of scattering 'islands' in the bulk. They describe the wires with a tight binding lattice model, in 
which some lattice sites are nonconducting (infinite potential). The conduction is then calculated numerically and 
averaged over a large number of configurations. Our model can be compared to their data for the case when there 
is no boundary roughness. In their fig. 8 our theory would predict a slightly higher average conductance (about 1 
times e'^/h higher) while the qualitative shape of the curve is about the same. It should be noted that a sample with 
these dimensions contains on the average only 5 (!) scatterers, and our extrapolation length zq is about 1.5 times the 
sample length dz- Therefore this sample is only in the crossover region between the ballistic regime and the regime 
of semi-ballistic transport. 

The data of Nikolic and MacKinnon also shows that the localization length for samples containing these strong 
scatterers is only a few times larger than the mean free path, which means that the range of validity of our analytic 
model, which presently does not contain localization corrections, is rather limited. Many experimentally accesible 
systems fall within this range though. 

Also, for weaker scatterers (with finite potentials) the localization length will be longer, extending the range of 
validity of our model to longer samples. 

Although our model assumes the scatterers are in the bulk of the wire, it can also be used to study qualitatively 
the effects of surface roughness, described as a high concentration of scatterers near the surfaces, with a vanishing 
scatterer density in the bulk. In this case there will be a strong subband bottom transparency effect when the second 
subband opens, as the scatterer density vanishes near the node of the second subband wavefunction. This gives rise 
to a peak in the transmission of the wire at the exact bottom of the second subband, as can be clearly seen in figures 
3 and 10 of |2l[. This peak is now well explained by our theory. 



VII. DISCUSSION 



We have found expressions for average transmission properties of two geometries of semi-ballistic devices, the 
waveguide-like geometry and the double barrier structure. Both geometries apply to optical systems and quantum- 
electronical systems. Partly such systems are already experimentally accessible, while other realizations are expected 
to become available soon. 

Steps in the conductance of electronic systems are predicted, together with remarkable drops shortly before a new 
subband opens. These drops have been described before in ref (surke). We notice that for systems containing pointlike 
scatterers, the behaviour of the conductance near these drops depends on the sign of the scattering length. 

Our expressions go beyond the second order Born approximation to the t-matrix of the individual scatterers (or 
random Gaussian potentials) and can therefore include scattering resonances. Divergencies are dealt with in a phys- 
ically meaningful way. It is shown that in a narrow waveguide the second order Born approximation is not always 
good, even when it is accurate for the same scatterers in wide threedimensional systems. 

Comparing our theory to numerical solutions of an Anderson model we find good agreement in the regime where 
localization effects can be neglected. 
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FIG. 1. The waveguide (with cross secion dxdy) is infinitely long but the disorder is only present in a finite section of length 
dz- The infinitely long 'pure' sections act as ideal quantum leads. 

FIG. 2. Conductivities of semi-ballistic devices with very low disorder, in the second order Born approximation. The solid 
line represents the conductivity of a disordered film of width d as it goes through multiples of the resonant width d*. The 
dashed line gives the conductivity of a ID constriction of a 2D electron gas as a function of the width d. 

FIG. 3. A double barrier quantum well, or its optical analogon, the Fabri Perot Interferometer, is modeled by a pair of 
delta function potential barriers in the x = and x = dx planes. Between these 'mirrors' impurities are present. Transport 
occurs in the x direction, tunneling through the barriers 

FIG. 4. Conductance of a long quantum wire containing a small density of strong attractive scatterers as a function of 
carrier energy E/V. The vertical lines indicate positions of the subband edges in the band structure of the lattice. 

FIG. 5. Conductance of a long quantum wire containing a small density of strong repulsive scatterers as a function of carrier 
energy. The vertical lines indicate positions of the subband edges in the band structure of the lattice. 

FIG. 6. Conductance of a long quantum wire with a high density of strong scatterers. At this point our simple model 
breaks down as localization corrections become important. The effective mean free path is in the order of one tenth of the wire 
length. 



FIG. 7. Extrapolation length zO and contact resistance Rc = zo/cr of a modestly disordered wire. 

FIG. 8. Conductance of short quantum wires with a large amount of attractive scatterers. The thick line represents the 
semi-ballistic theory, the thin line is an average over 10 numerical simulations. 

FIG. 9. Conductance of short quantum wires with a large amount of repulsive scatterers. The thick line represents the 
semi-ballistic theory, the thin line is an average over 10 numerical simulations. 

FIG. 10. Conductance of a 1000x11 quantum wire with 5% attractive scatterers, calculated using the full t-matrix (dashed 
curves) or the second order Born approximation(dotted curves). For a bare scattering strength of 0.03 the second order Born 
approximation is a rather good approximation to the full t-matrix, for u0=0.2 there is a large difference. 
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